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ABSTRACT 

I present the first results from a HST/NICMOS imaging study of the most metal-poor blue compact 
dwarf galaxy IZwl8. The near infrared color-magnitude diagram (CMD) is dominated by two popula- 
tions, one 10-20 Myr population of red super giants and one 0.5-5 Gyr population of asymptotic giant 
branch stars. Stars older than 1 Gyr are required to explain the observed CMD at the adopted distance 
of 12.6 Mpc, showing that IZwl8 is not a young galaxy. The results hold also if the distance to IZwl8 
is significantly larger. This rules out the possibility that IZwl8 is a truly young galaxy formed recently 
in the local universe. 

Subject headings: galaxies: dwarf, galaxies: compact, galaxies: stellar content, stars: AGB and 
post-AGB, infrared: galaxies, galaxies: individual (IZwl8) 



1. INTRODUCTION 

The blue compact galaxy (BCG) IZwl8 (also known as 
Markarian 116) was first described by Zwicky in 1966. In 
a seminal paper Searle and Sargent ( 1972| ) showed that 
IZwl8 had a very low oxygen abundance, and they con- 
cluded that the galaxy either must be young or form stars 
in short bur sts in tervened by long quiescent epochs (cf. 



an interesting idea that deserves to be tested. With its 
record low oxygen abundance, IZwl8 provides the critical 
test case. 

The chemical abundances of IZwl8 have been used to 



Searle et al. 1973). This has been followed by many stud- 



ies of the chemical p roper ties of IZwl8 and BCGs in gen- 
eral. Alloin et al. ( 1978 ) derived an oxygen abundance 
of 12+log(0/H)=7.2 (2% solar) for IZwl8 in agr eemen t 
with more recent values, e.g. Izotov and Thuan ( |1999|) . 
In spite of many searches for more metal-poor galaxies, 
IZwl8 remains the most metal-poor one (when nebular 
oxygen abundances are considered) and is often regarded 



as the prototypical BCG (see Kunth and Ostlin 2000 for 
a review). 

If IZwl8 was a genuinely young galaxy, that would im- 
ply that galaxies can form still at the present cosmic epoch 
and considerable effort has therefore been devoted to ex- 
amine the ages of BCGs. Although it is now clear that 
most BCGs are old (see Kunth and Ostlin |2000| ), IZwl8 
has remained among a small number of galaxies which 
show no clear signs for the presen ce of an old stellar popu- 
lation. Kunth and Sargent ( 1986j ) argued that IZwl8 could 
be young and that the low, but non-zero, metallicity could 
be due to self-pollution by massive stars. Recently an ex- 
tensive study of chemical ab undan ces of BCGs was pub- 
lished by Izotov and Thuan (199E). They found that for 
very low oxygen abundances (less than 5% solar) the N/O 
and C/O ratios stay constant and have very little scatter. 
Their interpretation is that N and C are produced entirely 
by massive (M > 10M Q ) stars and that this implies that 
these galaxies must be younger than 40 Myr. While other 



interpretations are possible (Kunth and Ostlin 2000), this 



infer both a high and a low age (Kunth and Ostlin 2000 ), 
and it is clear that chemistry alone has limited value 
in constraining the age and star formation history given 
the many uncertain physical parameters involved (winds, 
yields, infall, initial mass function, etc.). Stronger con- 
straints would be obtained by direct observations of the 
resolved stellar population. Hunter and Thronson ( 1995 ) 
used deep HST/WFPC2 imaging to investigate the stel- 
lar population of IZwl8, and found a population of young 
massive stars. No evidence for old stars was found, but 
neither any again st. Similar results were obtained by Du- 
four et al. ( 1996 ) from an independent WFPC2 data set. 
Lately, these two WFPC2 data sets hav e bee n carefully 
reanalyzed by Aloisi, Tosi and Greggio ( 1999 ), who dis- 
covered a population of red stars. Their results indicate 
that IZwl8 is at least .5 Gy r old. This was however dis- 
puted by Izotov et al. ( 1999 ) who claim that the distance 
to IZwl8 has been underestimated, and that hence the age 
has been overestimated. 

In this paper I report on the first results from a near 
infrared imaging study of IZwl8 using NICMOS on board 
the HST. The goal was to search for an old stellar popu- 
lation, by identifying luminous red giants and asymptotic 
giant branch (AGB) stars. The data is however also valu- 
able for studying the young massive stars, e.g. red super 
giants (RSGs) in the near infrared (NIR). A follow up 
program designed at investigating the presence of carbon 
stars by targeting molecular features were performed dur- 
ing Cycle 7-NICMOS. The results of this investigation will 
be presented in a forthcoming paper. 

2. OBSERVATIONS AND REDUCTIONS 
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HST/Nicmos imaging of IZwl8 



I Zw 18 was observed for 3 orbits during Cycle 7 with 
NICMOS in the F110W and F160W passbands (roughly 
corresponding to J and H) under GO program 7461. The 
instrument and its performance are described in Thomp- 
son et al. dl997D a nd th e NICMOS instrument handbook 



(MacKenty et al. 1997). The NIC2 camera was chosen 
to give good pixel sampling (0.075 "/pixel) and an ac- 
ceptable field of view (19" x 19"). All science exposures 
were obtained with the MULTIACCUM mode which em- 
ploys consecutive non-destructive readouts of the detector. 
This enables e.g. an efficient cosmic-ray removal. For the 
F110W filter, two MIF1024 sequences were acquired, sep- 
arated by an offset of 10 integer pixels in the Y direction. 
In addition, one MIF512 sequence was obtained with a 0.5 
pixel offset in both the X and Y directions with respect 
to the first MIF1024 sequence. The use of a non-integer 
pixel offset was adopted since the NIC 2 pixel scale slightly 
undersamples the point spread function (PSF) in F110W. 
For the F160W filter, two sets of images, consisting of one 
MIF512 and two MIF1024 sequences, were acquired at the 
same positions as the two MIF1024 sequences in F110W. 
The total exposure times are thus 2560 and 5120 seconds 
in the F110W and F160W filters respectively. 

Inspection of the pipeline processed images revealed that 
these suffered from the so called "pedestal" effect, a resid- 
ual bias not removed by the pipeline. To subtract the 
pedestal, the Pedestal Estimation and Quadrant Equaliza- 
tion Software developed by R.P. van der Marel was suc- 
cessfully utilized. After the pedestal had been removed, 
the images in each pass band were comb ined using the 
Drizzle package (Fruchter and Hook 1998). The resulting 
F160W image is shown in Fig. [[} 

Photometry of point sources was performed with the 
daophot package (Stetson et al. |1990| ). Unfortunately, the 
images contain few suitable stars for modeling the point 
spread function (PSF), since the bright stars are found in 
crowded regions. One bright source 4" south-west (SW) of 
the main body turned out to be spatially extended and is 
probably a background galaxy. To assure that the wings of 
the PSF were not contaminated by neighboring stars, pho- 
tometry was also performed using a synthetic PSF, which 
was constructed using the Tiny Tim software (version 4.4, 
Krist and Hook 1997| ). The resulting photometry using 
real and synthetic PSF stars were nearly identical, except 
for a systematic difference of the order of 0.05 magnitudes 
in the F160W filter; the use of synthetic PSF stars yielded 
redder stars. Comparison with a control sample consisting 
of the eight best isolated stars showed that the synthetic 
PSF photometry did not suffer from any systematic offset 
compared to the aperture photometry, whereas the use of a 
natural PSF showed a small systematic offset but smaller 
scatter than the synthetic PSF photometry. Thus, the 
natural PSF photometry is probably more accurate but 
subject to a small systematic offset. No systematic dif- 
ference was seen in the F110W data. The natural PSF 
photometry was finally adopted after correcting for the 
small offset in F160W. After the photometric reduction, 
stars that resulted too sharp or too extended in at least 
one filter were rejected. The photometries in the two filters 
were then combined by considering a matching radius of 
0.5 pixels in X and Y for the stellar centers in the F110W 
and F160W images. All objects were carefully analyzed 
by eye and dubious sources rejected. 



The DAOPHOT magnitudes refer to an aperture with 
radius 0. 1875". The aperture corrections in going from this 
to an 0.5" aperture were determined using synthetic PSF 
and amount to 0.24 and 0.48 magnitudes for F110W and 
F160W respectively. To obtain total magnitudes an aper- 
ture correction of —0.156 magnitudes for both filters were 
added to the the 0.5" aperture magnitudes, and finally, the 
zeropoints of ZPpxiow = 22.5436 and ZPpxmw = 21.826 
were added in order to have the final calibrated magni- 
tudes in the HST Vegamag system (see NICMOS WWW 
pages). I will use the designations mno and TOi6o to ex- 
press the magnitudes in F110W and F160W respectively, 
and Fl 10W - F160W for the color. The resulting color- 
magnitude diagrams (CMDs) and the magnitude-error re- 
lations are shown in Fig. [| The detected stars span a 
similar range in color and brightness as those detected 
with HST/NICMOS in VIIZw403 (Schulte-Ladbeck et al. 
1999| ). 

The photometric completeness was investigated by using 
the addstar task. Artificial stars were added in half mag- 
nitude bins at random locations in two separate regions: 
"SE" covering the mildly crowded south-east region; and 
"NW" which is the severely crowded north-west region. 
Then, repeated tests were made of how many test stars 
within each half magnitude bin that were successfully re- 
covered by daophot. More than 100 stars were added for 
each bin. The 75 and 50 % completeness levels are sum- 
marized in Table 1. In the NW region there are still bright 
unresolved sources left after the repeated daophot runs. 

3. THE DISTANCE OF IZW18 

Usually a distance of 10 Mpc is assumed for IZwl8, 
based on a radial velocity of ~750 km/s and a Hubble 
consta nt of Hq = 75km/s/Mpc. Recently Izotov et al. 
have claimed that the distance to IZwl8 has been 
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severely underestimated, and if Virgo-centric flow is taken 
into account a distance of 15-20 Mpc is obtained. They 
present also other circumstantial evidence (the presence of 
Wolf-Rayet stars and the ionization state of IZwl8) for a 
distance D ~ 20 Mpc. It is true that Virgo-centric mo- 
tion will slightly increase th e dist ance, but the distance 
proposed by Izotov et al. ( 1999 ) refer to old values of 
the distance of the Virgo cluster (21.7 Mpc). A lot of 
recent progress has been m ade on the distance to Virgo 



e.g. Macri et al. 1999 ) showing that its distance is 



Dvirgo ~ 16 Mpc. The Virgo-centric flow models express 
the distance of a galaxy as x, relative to the Virgo clus- 
ter, i.e. D = x x Dvirgo ■ Applying the linear model by 



Schcchtcr (198C) gives x = 0.7 while the non-linear mod- 
els by Kraan-Korteweg ( |1986| ) give x — 0.8, correspond- 
ing to 11 and 13 Mpc respectively. Hereafter I will adopt 
D = 12.6 Mpc for IZwl8, which corresponds to a distance 
modulus of (to — M) = 30.5. The distance should not be 
much larger (maximum 14.5 Mpc given the Virgocentric 
flow models and the uncertainty in Dvirgo) than this, but 
may be as low as ~ 10 Mpc. The distances suggested by 
Izotov et al. ( 1999 ) would require IZwl8 to be more dis- 
tant than the Virgo cluster, in spite of a much lower radial 
velocity. 



4. ANALYSIS AND DISCUSSION 
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Internal extinction has only a moderate influ- 
ence on the F110W-F160W colors and in general 
£(F110W - F160W) = 0.5 x E(B - V), which has been 
estimated with the Synphot software. All optical investi- 
gations have found a low internal extinction in IZwl8. I 
have adopted E{B — V) = 0.05 throughout the galaxy, and 
this correction has been applied to the data-points in Fig. 



The obtained color magnitude diagram (CMD) clearly 
presents two distinct features. One population of red lu- 
minous stars and one with fainter (mno > 23) red stars. 
This suggests that stars with at least two different ages 
are present. To interpret the observed CMD I have mainly 
used the Padua stella r models for a metall icity Z = O.O2Z 
(Fagotto et al. |1994| , Bertelli et al. |1994[ ), i.e. identical to 
the ga s metallicity of IZwl8. In general, the Bertelli et al. 
( 1994 ) isochrones are more useful than the Fagotto et al. 
(1994) tracks since the former contains later evolutionary 
stages, as the thermally pulsating AGB phase, which is 
very important for the present study. 

Before interpreting the CMD, the effective temperature 
(T G ff) and luminosity (L) of the stellar models have to be 
related to colors and absolute magnitudes in the observed 
system. To translate T e g into F110W-F160W colors, and 
L into Mno for a given surfac e gravity, the BCP models 
by Origlia and Leitherer (2000) were used for a metallicity 
of [Fe/H] = -2. 

In figure || the Bertelli et al. (1994) isochrones are 
over-plotted on the CMD. For a distance modulus of 



(m — M) = 30.5 it is obvious that stars of ages 11-20 
Myr are required to reproduce the luminous upper branch 
in the CMD (m 110 < 23.2, F110W - F160W < 1). These 
ages equals the lifetime of stars with initial masses of 20 to 
13 Mq. The upper branch is dominated with stars « 14 
Myr old when the current star formation event appears 
to have peaked. For comparison, distances of 10 and 15 
Mpc would give ages of 18 and 1 2 My r respectively. Com- 
parison with the Fagotto et al. ( 1994| ) stellar evolutionary 
tracks yield the same conclusion that the most luminous 
red stars are the evolved descendants of stars with initial 
masses 15-20 Mq, and the same results are obtained when 
using "Geneva" tracks for the same metallicity (Scha erer 
1999, private communication; cf. de Mello et al. 199S). 



For low mass stars (initial mass M <J 5M Q , age > 100 
Myr) only the isochrones (Bertelli et al. 1994) have 
been considered, since the tracks do not include the ther- 
mally pulsating AGB phase. The most interesting re- 
sult is that to reproduce the faint red part of the CMD 
(mno > 23.5,F110VF - F160W > 1), an age of 1 Gyr 
is required at the distance of 12.6 Mpc. Even if the dis- 
tance would be as high as 15 Mpc, stars older than 0.5 
Gyr would be required, and D — 10 Mpc would lead to 
a much higher age than 1 Gyr. Not even at 20 Mpc can 
the faint red stars be made younger than 300 Myr. For 
D = 12.6 Mpc the red faint population can be nicely ex- 
plained as the upper part of AGB stars with a range of 
ages from 0.3 to 5 Gyr. This corresponds to the lifetime of 
stars with initial masses 2.5 to 1.1 Mq. A single 0.5 Gyr 
old AGB population does not suffice since that would not 
produce enough faint (mno ~ 24) red stars, only when 
extending to ages in excess of 1 Gyr can the CMD of faint 
stars be satisfactory explained. This indicates that IZwl8 



has experienced more or less continuous star formation for 
at least 1 Gyr prior to the present burst. 

The distribution of stars in the CMD is broadened by 
the observational error function, and it has to be assured 
that it cannot be reproduced by young stars plus pho- 
tometric errors. If one assumes that no star is older 
than 100 Myr, what distribution would one expect? An 
isochrone of that age has a color F110W — F160W = 0.79 
at mno = 24. The one sigma photometric errors at mno = 
24, F110W - F160W w lare (Tfiiqw-fwow = 0.125, but 
let's be very conservative and adopt (Tfuow-fwow = 0-2. 
The one sigma range around F110W — F160W = 0.79, 
mno ~ 24 then contain 14 stars. For a normal distri- 
bution, 68% of the stars should be contained within this 
range, and one would thus expect to find only ~ 3 stars 
(i.e. 16%) with FllQW - F160W > 1.0, wheras one finds 
more than a dozen. Thus the observed distribution of red 
stars cannot be reproduced by 0.1 Gyr old stars plus pho- 
tometric errors. Even if rejecting the reddest stars (> 1.5), 
an age of 1 Gyr is required to be consistent with the er- 
ror distribution. Since this was a conservative estimate, 
we can conclude that IZwl8 has an age of at least 1 Gyr. 
This result on the near-infrared CMD analysis is confirmed 
when using the synt hetic CMD technique (cf. Tosi et al. 
(1991); Aloisi et al. 1999): stars older than 1 Gyr are re- 
quired in order to reproduce the observed CMD, the best 
fit being obtained by including stars as old as 5 Gyr (Ostlin 
2000, in preparation). 

An analysis of the spatial distribution of stars of dif- 
ferent age show that young red super giants occupy pref- 
erentially the NW region, while the old AGB stars are 
more evenly distributed. Further constraints on the age 
of the stellar population may be obtained from the color 
of the underlying galaxy after the found stars have been 
subtracted. For a single stellar population, the integrated 
F110W-F160W color increases with age, from F110W — 
F160W = -0.1 at 10 Myr up to F110W - F160W = 0.6 
at 125 Myr when the AGB stars start to dominate, and 
stays roughly constant thereafter. This has been calcu- 
lated using the PEGASE spectral synthesis code (Fioc & 
Rocca-Volmerange 1999) and the NICMOS F110W and 
F160W response curves. The integrated colors in 15 x 15 
pixels boxes are typically F110W - F160W = 0.3 to 0.6 
in the NW and F110W - F160W > 0.6 in the SE re- 
gion. This indicates that in the NW region, the light is 
dominated by an unresolved population of stars with ages 
20 to 100 Myr, while in the SW region, the underlying 
population is dominated by stars older than 125 Myr. 

The result of an age of 1 Gyr or larger for the underly- 
ing population in IZwl8 is confirmed also by deep optical/ 
near-IR surface photometry (Ostlin 2000, in preparation; 
Kunth & Ostlin |2000| ) a nd is in good agreement with the 
results of Aloisi et al. ([19990 for the optical CMDs. We 
can thus conclude that all these different studies on inde- 
pendent sets of data in various spectral regions give the 
same answer: that IZwl8 is an old galaxy. 

5. SUMMARY AND CONCLUSIONS 

I have analyzed the near infrared color-magnitude di- 
agram from HST/NICMOS images of IZwl8, the most 
metal-poor galaxy known. Assuming a distance of 12.6 
Mpc, I find support for the presence of AGB stars at least 



HST/Nicmos imaging of IZwl8 



1 Gyr old. The results presented in this letter hold qual- 
itatively even if the distance to IZwl8 is significantly dif- 
ferent. The latest star formation event appears to have 
peaked 14 Myr a go. M y results are in good agreement 
with Aloisi et al. (1999) who from a careful reanalysis of 



HST/WFPC2 archival data concluded that IZwl8 must 
be at least 0.5 Gyr old at their adopted distance of 10 
Mpc. Thus, there is no need for galaxies with oxygen 
abundance less than 5 % of the sola r valu e to be young as 
was claimed by Izotov and Thuan (1999). This closes the 
question whether IZwl8 is a genuinely young galaxy, still 
in the process of formation: the answer is - No! 
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HST/Nicmos imaging of IZwl8 




Fig. 1 — The central 14.5" X 14.5" of the drizzled NICMOS F160W image of IZwl8. The north and cast directions are indicated at the 
lower right. The length of the arrow is 2". 
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Fig. 2. — The left panels show the observed color-magnitude diagrams, and the right panels show the observed magnitudes versus 
photometric uncertainties for F110W (top) and F160W (bottom). 
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HST/Nicmos imaging of IZwl8 




F110W-F160W 



Fig. 3. — Observed color— magnitude diagram (corrected for extinction) with isochrones from Bertelli et al. (1994) over-plotted. The 
numbers correspond to the logarithm of the age for each isochrone. The right hand axis gives the absolute F110W magnitude for a distance 
of 12.6 Mpc. 
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Table 1 

Completeness limits in magnitudes, for the two regions discussed in the text. 



